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ABSTRACT 


A two-dimensional mean field dynamo model is solved 
where magnetic helicity conservation is fully included. 
The model has a negative radial velocity gradient giving 
rise to equatorward migration of magnetic activity pat¬ 
terns. In addition the model develops longitudinal vari¬ 
ability with activity patches travelling in longitude. These 
patches may be associated with active longitudes. 

Key words: Magnetohydrodynamics (MHD) - turbu¬ 
lence - Sun: active longitudes. 


1. INTRODUCTION 


Active regions are complexes of magnetic activity out of 
which sunspots, flares, coronal mass ejections, and sev¬ 
eral other phenomena emerge with some preference over 
other regions. These regions tend to be bipolar, i.e. they 
come in pairs of opposite polarity and are roughly aligned 
with the east-west direction. 

There is some controversy as to whether or not active 
regions appear preferentially along so-called active lon¬ 
gitudes and what their long term stability properties are 
(e.g., Bai 2003, Berdyugina & Usoskin 2003, Pelt et 
al. 2005). Some degree of recurrence of sunspots has 
frequently been reported over the years (Vitinskij 1969, 
Bogart 1982, Bai 1987, 1988), but only now these ideas 
are becoming more quantitative. 

As the recent analysis of Usoskin et al. (2005) has shown, 
active longitudes have characteristic angular velocities 
that depend on the phase of the cycle which, in turn, 
determines the typical latitude of their occurrence. The 
analysis of solar magnetograms by Benevolenskaya et al. 
(1999) showed already that at the beginning of each cy¬ 
cle, when most of the activity occurs at about ±30° lati¬ 
tude, the rotation rate of the active longitudes is Cl/2n ss 
446 nHz, while at the end of each cycle, when the typical 


latitude is only ±4° latitude, the rotation rate of the ac¬ 
tive longitudes is Q/2n ss 462 nHz. The recent work of 
Usoskin et al. (2005) demonstrates that the active regions 
can also be detected in sunspot data. Unlike the analysis 
of Benevolenskaya et al. (1999), who isolated two dif¬ 
ferent active longitudes at the beginning and the end of 
the cycle, Usoskin et al. (2005) determined a continuous 
latitudinal dependence of these active longitudes on the 
phase of the cycle. The success of their analysis lies in 
the way they calculated time-dependent reference values. 
In particular, they find that only 10% of the spots partic¬ 
ipate in this nonaxisymmetric effect. Thus, the effect is 
real, but weak, although with a well determined strength. 

Given that active longitudes suggest the presence of 
well preserved activity patches in an otherwise turbulent 
medium, one has to look for a quantity that has the capa¬ 
bility to be long-lived. An obvious candidate for anything 
long lived in hydromagnetic turbulence at large magnetic 
Reynolds numbers is the magnetic helicity. In the ab¬ 
sence of magnetic helicity fluxes, magnetic helicity is 
nearly perfectly conserved in resistive magnetohydrody¬ 
namics at large magnetic Reynolds numbers. Even in the 
presence of magnetic helicity fluxes, magnetic helicity is 
only transferred from one place to another, but it is not 
lost. It is therefore plausible that there might be a con¬ 
nection between the nearly perfect magnetic helicity con¬ 
servation and the long-lived features on the sun such as 
active longitudes. In the simplest case, magnetic helicity 
could be considered frozen into the plasma so that local 
patches of enhanced magnetic helicity would just propa¬ 
gate with the ambient velocity of the gas. This could in 
principle be a very simplistic picture of active longitudes 
that might explain their long life times, but not how they 
came into existence. 

An alternative that is traditionally discussed in connec¬ 
tion with active longitudes is the idea that there are a 
number of different axisymmetric and nonaxisymmet¬ 
ric dynamo modes present in the sun that are mixed in 
the right proportions such that their superposition corre¬ 
sponds to the observed field configuration (Radler et al. 
1990, Moss 1999, 2004, Moss & Brooke 2000, Bigazzi 
& Ruzmaikin 2004, Berdyugina et al. 2006). An obvious 


problem is that these modes are solutions of the linearized 
problem and that their superposition does not constitute a 
solution to the nonlinear problem. At first instance only 
one of the modes gets selected, so the final solution is 
sill mainly either axisymmetric or nonaxisymmetric, but 
not easily anything in between, as originally anticipated 
(Radler et al. 1990, Moss et al. 1995). 

Magnetic helicity conservation only applies to the total 
held, i.e. the sum of small scale and large scale fields. 
The large scale held produced by a mean held (large 
scale) dynamo does by itself not conserve magnetic he¬ 
licity, because the a effect leads to a transfer of magnetic 
helicity from smaller to larger scales (Pouquet et al. 1976, 
Ji 1999). It is therefore necessary to include also the con¬ 
tribution from the small scale magnetic helicity, which 
enters into the mean held description through a magnetic 
contribution to the a effect. This approach is now fairly 
well developed and is important in reproducing the slow 
saturation (Field & Blackman 2002, Blackman & Bran¬ 
denburg 2002, Subramanian 2002) found in simulations 
(Brandenburg 2001). For a review of these recent devel¬ 
opments see Brandenburg & Subramanian (2005a). 

We begin by explaining this approach in more detail. In 
order to focus on the essentials, we consider only a min- 
imalistic model. The basic dynamo wave can already be 
described in a one-dimensional model with only latitudi¬ 
nal extent. In the present context we still need the longi¬ 
tudinal extent, but we ignore the variation of the magnetic 
held with depth. Furthermore, in order to study the basic 
effect introduced by the longitudinal extent and by mag¬ 
netic helicity conservation, it suffices to restrict oneself 
to a cartesian model. Since there is no radial extent, there 
are also no radial boundaries, and hence no magnetic he¬ 
licity flux out of them. This might be an important limi¬ 
tation to keep in mind. 


2. THE MODEL 


The dynamo equation together with the dynamical 
quenching model is (Field & Blackman 2002, Blackman 
& Brandenburg 2002, Subramanian 2002) 


scale, D/D t = d/dt + U • V is the advective deriva¬ 
tive with respect to the mean How U , and 

£ = q(aB - p t p, 0 J) (4) 

is the turbulent electromotive force (assuming local 
isotropy), with an ad hoc added overall quenching of the 
electromotive force via the factor 

q=l/(l + q B B 2 /B 2 eq ). (5) 

Such an overall quenching was applied also in the recent 
dynamically quenched mean held of Brandenburg & Sub¬ 
ramanian (2005b). 

Under the assumption of locally isotropic turbulence we 
have a = qk + oh. where «k = — • u is the 

kinetic a effect, where ui = V x u is the vorticity, 
a m = ^rj-b/pois the magnetic contribution to the a 

effect, and ry = \tu 2 is turbulent magnetic diffusiv- 
ity. In deriving equation □ we have used the defini¬ 
tions £?g q = popu 2 and rj t = \tu 2 , so that ^r/(pop) = 
/i/q/r/ t (see Blackman & Brandenburg 2002). 

Equation 0 can now be solved together with equa¬ 
tion m to determine the effect of magnetic helicity evo¬ 
lution on the solar dynamo. These equations have been 
solved in recent years both in context of the galactic dy¬ 
namo (Kleeorin et al. 2000, 2002, Shukurov et al. 2005) 
and the solar dynamo (Kleeorin et al. 2003, Zhang et al. 
2006). The dynamical quenching model has also been 
tested against direct simulations (Blackman & Branden¬ 
burg 2002, Brandenburg & Subramanian 2005b). 

We emphasize however that, although B depends only 
on x and y, U depends also on z in a prescribed fashion, 
such as to model the radial differential rotation. To en¬ 
sure that the held stays numerically divergence free, it is 
advantageous to write equation 0 in terms of the mean 
magnetic vector potential A, so B = V x A and the ra¬ 
dial velocity gradient is explicitly added to the model by 
working with the velocity in the form 

U = U { °\z) +u(x,y,t)-, U {0) = (53,0,0). (6) 
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where B is the mean magnetic held, J = V x B / p,Q is 
the mean current density, po is the vacuum permeability, 
77 is the microscopic magnetic Spitzer diffusivity, ry is the 
turbulent magnetic diffusivity. 


Rm = Vt/v 


(3) 


We apply the equations at the reference height 2 = 0, so 
U < ' 0 ' 1 itself is zero and does not enter the equations; only 

a single component of its derivative matrix, U^\ = S, 
enters. Thus, the full set of equations is 
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is here used as the dehnition of the magnetic Reynolds 
number, fcf is the wavenumber of the energy-carrying 
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Figure 1. Time dependence of the rms mean magnetic 
field and the rms mean velocity, showing the emergence 
of oscillatory and later more irregular solutions. 

We model the magnetic field only in one hemisphere and 
restrict ourselves to dipolar parity by demanding 

A x ,y = Ay = A Z y = 0 (on y = 0), (11) 

corresponding to a normal-field condition (n x B = 0), 
modeling the equator, and 

A x = Ay. y = A z = 0 (on y = 2tt), (12) 

corresponding to a perfect conductor boundary condition 
(n ■ B = 0), in an attempt to capture some of the behav¬ 
ior on the pole. The simulations have been carried out 
using the PENCIL CODE 1 which is a high-order finite- 
difference code (sixth order in space and third order in 
time) for solving the compressible hydromagnetic equa¬ 
tions. The code comes with a special mean field module 
for the dynamical a quenching equations. 




3. RESULTS 


In the following we discuss in detail one particular model 
using the following set of parameters: S = —30 (corre¬ 
sponding to negative radial differential rotation), qk = 
0.5, t't = rjt = 1, t] = 0.05 (so R m = 20), k{ = 5, and 
an overall quenching factor with qn = 0.1. 

In Fig. [0 we show the evolution of the rms values of the 
mean magnetic field, (B ) 1 / 2 and of the mean velocity 

field, (U ) l,/2 . Initially, the magnetic field grows expo¬ 
nentially, but then it saturates in an oscillatory fashion. 
The mean field also drives a mean flow. The long time 
behavior is more complicated, as is shown in the inset: 
the velocity increases slightly, but then the field becomes 

1 http://www.nordita.dk/software/pencil-code 


Figure 2. Time dependence of various integral quanti¬ 
ties, in the quiescent phase after rjtkito = 2700. The 
vertical bars indicate the times for which synthetic ‘mag¬ 
netograms’ and plots of other quantities are given below. 


more irregular. There are also stages where the veloc¬ 
ity exceeds the magnetic field and limits the time depen¬ 
dence. 

In Fig. |3 we show the evolution of the rms values of ve¬ 
locity and magnetic field, as well as the maximum veloc¬ 
ity, together with the mean magnetic and current helici- 
ties, (A ■ B) and (J ■ B) during a short time interval dur¬ 
ing a quiescent stage. The vertical bars mark the times 
for which snapshots of various quantities are shown in 
Figs [30 that will be discussed next. 









































Figure 3. Plot of B z (x 1 y) at different times, in regular 
intervals ofr] t k\At = 0.5 after t = to- 


We begin with Fig.0 where we show images of the line 
of sight component of the magnetic field, B z . One clearly 
sees magnetic activity patches moving both equatorward 
(downward in the plot) as well as in the prograde direc¬ 
tion (to the right in the plot). The fact that these patches 
propagate at all is interesting. It is not related to a locally 
enhanced mean flow in the prograde direction: the flow 
is actually in the retrograde direction. The propagation 
must therefore be related to a proper three-dimensional 
dynamo mode travelling to the right, just like in the non- 
axisymmetric dynamo models mentioned in the begin¬ 
ning. 

Next, we consider images of current helicity, J ■ B, in 


Figure 4. Plot of J ■ B at different times, in regular 
intervals ofrjtkiAt = 0.5 after t = to- 


Fig.0 which clearly shows enhanced positive current he¬ 
licity in the regions where also the line-of-sight magnetic 
field is strong. The current helicity has the same sign in 
the magnetic patches with positive and negative net flux. 

Finally, we consider images of the magnetic a effect. The 
connection between q:m and field strength or current he¬ 
licity in not very clear. Indeed, the two are only loosely 
connected with each other. The primary reason for hav¬ 
ing the magnetic a term is to provide a nonlinear feed¬ 
back between the produced large scale helicity and the a 
effect that helps to maintain total magnetic helicity con¬ 
servation. 












Figure 5. Plot of ayi{x, y) at different times, in regular 
intervals ofrjtkiAt = 0.5 after t = to- 


4. CONCLUSIONS 


The present investigations can only be regarded as pre¬ 
liminary, because we have not been able to explore the 
big parameter space given by the large number of un¬ 
knowns. Most worrisome is the restriction to only small 
values of the magnetic Reynolds number. At the moment 
we are experiencing problems when we try larger val¬ 
ues, suggesting that the problem may not we well-posed 
for larger values of R m and may require modifications. 
However, this is still surprising, because similar problems 
have not been encountered in simpler problems, where 
however no mean advection or shear was taken into ac¬ 


count (e.g. Brandenburg & Subramanian 2005b). 

One surprising aspect emerging from the present simula¬ 
tions is the fact that magnetic patterns do not move with 
the local gas velocity. Instead, the field propagation is 
governed by the nonaxisymmetric dynamics of the dy¬ 
namo modes. Thus, these simulations would not support 
the notion of field line anchoring. This pictures is oc¬ 
casionally used in connection with sunspot proper mo¬ 
tions. Long before the internal angular velocity was de¬ 
termined via helioseismology, it was known that sunspots 
rotate faster than the surface plasma (Howard et al. 1984). 
Moreover, young sunspots rotate faster than old sunspots 
(Tuominen & Virtanen 1988, Pulkkinen & Tuominen 
1998). A common interpretation is that young sunspots 
are still anchored at a greater depth than older ones, and 
that therefore the internal angular velocity must increase 
with depth (see also Brandenburg 2005). This provided 
also the basis for the classical mean field dynamo theory 
of the solar cycle according to which the radial angular 
velocity gradient has to be negative (Parker 1987). 

Given the lack of agreement between the speed of mag¬ 
netic patches and the local flow speed, one it led to be¬ 
lieve that magnetic structures can therefore not be used 
as tracers of the local flow speed. Alternatively, it is also 
possible that the tracer properties of the advected field 
are fully displayed only in three dimensions, or at larger 
magnetic Reynolds numbers. However, the present re¬ 
sults should be interpreted with care, because the present 
calculations have only been possible in a rather limited 
parameter range and for rather small magnetic Reynolds 
numbers. It would be important to test the notion of field 
line anchoring in direct simulations of the original non- 
averaged equations, i.e. in the presence of developed tur¬ 
bulence. 
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